Astronomy & Astrophysics manuscript no. ha6'lc © ESO 2013 

February 24, 2013 



Integral field spectroscopy based Ha sizes of local luminous and 

ultraluminous infrared galaxies 

A direct comparison with higli-z massive star forming galaxies 

S. Arribas^, L. Colina^, A. Alonso-Herrero^'^, F.F. Rosales-Ortega^'^, A. Monreal-Ibero"* , M. Garcia-Marin^ , S. 

Garcia-Burillo^, J. Rodriguez-Zaurin' 

04 ' ' Centre de Astrobiologia, Departamento de Astrofisica, CSIC-INTA, Ctra. de Ajalvir km. 4, 28850 -Torrejon de Ardoz, Madrid, 

Spain. 

^D . e-mail: arribas@cab.inta-csic.es 

CN ■ - Institute de Fi'sica de Cantabria, CSIC-UC, Avenida de los Castros S/N, 39005-Santander, Spain 

» I ' ^ Departamento de Fisica Teorica, Universidad Autonoma de Madrid, 28049 Madrid, Spain. 

C^ , "^ Instituto de Astrofisica de Andalucia. CSIC. Glorieta de la Astronomia, s/n, 18008, Granada, Spain. 

' Physikalisches Institut, Universitat zu Koln, Ziilpicher Strasse 77, 50937 Koln, Germany 

* Observatorio Astronomico Nacional (OAN)-Observatorio de Madrid, Alfonso XII, 3, 28014-Madrid, Spain 

(^ ■ ' Space Telescope Science Institute, 3700 San Martin Drive, MD2I2I8, Baltimore, USA 

m ' 

ABSTRACT 

^^ ' Aims. We study the analogy between local luminous and ultraluminous infrared galaxies (U/LIRGs) and high-z massive star forming 

Vi^ I galaxies (SFGs) by comparing their basic Ka structural characteristics, such as size and luminosity surface density, in an homoge- 
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neous way (i.e. same tracer, size definition, and similar physical scales). 

Methods. We use integral field spectroscopy (IPS) based Ha emission maps for a representative sample of 54 local U/LIRGs (66 
galaxies) observed with INTEGRAL/WHT and VIMOS/VLT. From this initial sample, we select 26 objects with similar Ha lumi- 
Q I nosities (L(Hq')) to those of massive (i.e. M^. ~ 10'" Mq or larger) SFGs at z ~ 2, and observed on similar physical scales. We then 

^ , directly compare the sizes, and luminosity (and SFR) surface densities of these local and high-z samples. 

j/3 ■ Results. The size of the Ha emitting region in the local U/LIRGs that we study has a wide range of values, with ri/2(Ha) from 0.2 

C^ ' kpc to 7 kpc. However, about two-thirds of local U/LIRGs with L,> > 10" "* Lq have compact Ha emission (i.e. ri/2 < 2 kpc). The 

comparison sample of local U/LIRGs also contains a larger fraction (59%) of objects with compact Ho- emission than the high-z sam- 
ple (25%). This gives further support to the idea that for this luminosity range the size of the star forming region is a distinctive factor 
when comparing local and distant galaxies of similar SFRs. However, when using Ha as a tracer for both local and high-z samples, 
the differences are smaller than those reported using a variety of other tracers. In the L(Ha) - L(Ha) surface density (E/y„) plane, 
most of the local U/LIRGs and high-z SF galaxies follow the same trend (i.e. higher luminosity for higher surface density) and cover 
a similar range, except for about 20-40 % of local U/LIRGs, which have a higher E;/„ by a factor of about 10. This is considerably 
smaller than the factors of 1000 or more reported in similar planes (i.e. L(TIR) versus "Ltir)- Despite of the higher fraction of galaxies 
with compact Ha emission, a sizable group (about one-third) of local U/LIRGs are large (i.e. ri/2 > 2 kpc). These are systems that 
rf) ' show evidence of pre-coalescence merger activity and are indistinguishable from the massive high-z SFGs galaxies in terms of their 

C^ , Ha sizes, and luminosity and SFR surface densities. 

Key words, galaxies - sizes - luminous infrared galaxies - integral field spectroscopy 



'k>< 1. Introduction (i.e. L,> ~ IO'^Lq). Although U/LIRGs are rare locally, stud- 
'V^ ' ies with the Spitzer Space Telescope have shown that they are 
^ ; Luminous and ultraluminous infrared galaxies (LIRGs: L,, s ^^^^j^ ^^^^ numerous at high-z and account for an increasingly 
- - - L[8 - 1000//m] - 10 Lq; ULIRGs: L,, > 10 Lq) are j j. fj.^^,^^^^^ ^f (^e total star formation density fraction (e.g. 
beheved to have an important role in our understanding of ^^^.^ ^^^^ j^^jf ^^ ^j^^ ^^^^j gp ^^^^. ^^ ^^2, Perez-Gonzalez 
galaxy evolution. They are systems of intense star foririation ^^ ^j^ 2005; although see Rodighiero et al. 201 1). The so-called 
(SF), whose spectral energy distributions (SEDs) are domi- sub-millimeter galaxies (SMGs; e.g. Small et al. 1997) are com- 
nated by dust thermal emission arising from the reprocessing ^^^^ ^-^^^^ ^^ ^ ^^^^ luminous counterpart of local U/LIRGs. 
of UV photons produced by young massive stars and/or ac- 
tive galactic nucleus (AGN) heating. The fraction of these sys- However, although many high-z galaxies meet the above 
tems with disturbed morphologies increases with luminosity, (luminosity) definition, the analogy between local and distant 
with most ULIRGs showing evidence of recent or on-going LIRGs and ULIRGs is under discussion. Several authors have 
merger events. By comparison, LIRGs seem to be a more het- reported that the SEDs of high-z U/LIRGs are similar to lower- 
erogeneous group with in many cases properties similar to those luminosity local systems suggesting that they are not high-z 
of isolated star-forming spirals, especially at low luminosities analogs but instead scaled-up versions of lower-luminosity lo- 
cal U/LIRGs (e.g. Pope et al. 2006, Papovich et al. 2007, Takagi 
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It has been suggested that the observed differences between 
local and distant U/LIRGs (of similar luminosity) could be due 
to a difference in metallicity and/or physical size (e.g. Rigby et 
al. 2008; Farrah et al. 2008). On the one hand, Engelbracht et 
al. (2008) find that the differences between the SEDs of local 
and z ~ 2 ULIRGs are qualitatively consistent with a difference 
in metallicity of a factor 1.5-2. On the other hand, several au- 
thors find that local U/LIRGs have a more compact structure 
than high-z populations (e.g. lono et al. 2009, Rujopakarn et al. 
201 1, and references therein). Elbaz et al. (201 1) also linked the 
differences in SEDs to the compactness of the star-forming re- 
gion. 

Using a variety of tracers, Rujopakarn et al (201 1) report that 
local U/LIRGs have smaller SF region sizes by up to factors of 
50 (i.e. more than 3 orders of magnitude in luminosity surface 
density) than other local and high-z SFGs. Therefore, they sug- 
gest that local ULIRGs and LIRGs belong to a rare population 
driven by a unique process. These results have brought a lot of 
attention to the determination of U/LIRG sizes (especially for 
local samples). 

Despite size being a fundamental property, its determination 
is often uncertain as it is affected by a number of observational 
and instrumental factors (e.g. reddening, resolution, etc). In ad- 
dition, different spectral features trace different galaxy compo- 
nents and physical mechanisms. This complicates the compari- 
son between local and distant objects, which are often selected 
and observed at different (rest-frame) wavelengths. Differences 
in methodology also add uncertainty to this comparison. 

This paper compares some structural properties of the Ho- 
emitting region in local U/LIRGs, such as sizes and luminos- 
ity surface densities, with those of high-z populations. We use 
integral field spectroscopy (IFS) based Ha emission maps ob- 
tained from our INTEGRAL/WHT and VIMOS/VLT observa- 
tions to derive these properties for a representative sample of lo- 
cal U/LIRGs (e.g. Colina et al., 2005; Garcia-Marin et al, 2009a; 
Rodriguez -Zauiin et al, 201 1, and references therein). We com- 
pare our results with those for high-redshift massive SFGs sam- 
ples observed with near-infrared IFS on similar linear scales, 
using the same tracer (e.g. Forster-Schreiber et al, 2009, 2011; 
Wright et al., 2009; Law et al., 2009; Wisniosky et al. 2011). 
Therefore, we follow similar methods for the local and distant 
samples to minimize the uncertainties in the relative compari- 
son. As we base our size measurements on reddened Ha maps, 
they refer to the extension of the "unobscured" Ho'-emitting re- 
gion. 

The paper is structured as follows. In Section 2, we briefly 
describe the IFS sample and the data used for the analysis. 
In Section 3, we first derive SF region sizes from Ha for the 
whole sample of U/LIRGs observed with INTEGRAL/WHT and 
VIMOS/VLT. From this group, we then select the subsample that 
allows us to make a homogeneous comparison with the high-z 
samples. The sizes and luminosity surface densities of distant 
and local U/LIRGs samples are confronted in Sec. 4. Finally our 
main conclusions are summarized in Section 5. Throughout the 
paper, we consider Hq = 70 kms 'Mpc ', Qa = 0.7, Qm - 0.3. 

2. The sample and observations 

The sample for the present study is formed by local LIRGs and 
ULIRGs for which we have obtained optical IFS data with the 
VIMOS-IFU (LeFevre et al. 2003) and the INTEGRAL (AiTibas 
et al. 1998) instruments. The objects were selected to form a 
representative sample of U/LIRGs covering all types of mor- 
phologies (isolated spirals, interacting pairs, mergers), nuclear 



excitation (HII, Seyfert, and LINER), and sampling the LIRG 
and ULIRG infrared luminosity range. The sample is incomplete 
in either flux or distance, since complete samples covering the 
U/LIRG luminosity range would be so numerous that their IFS 
observations would require a prohibitively large amount of time. 

The INTEGRAL sub-sample is that presented in Garcia- 
Marin et al. (2009a). These are northern objects selected from 
the IR-bright samples of Sanders et al. (1988), Melnick & 
Mirabel (1990), Leech et al. (1994), Kim et al. (1995), Lawrence 
et al. (1999), and Clements et al. (1996). For the present analy- 
sis we exclude IRAS F09427H-1929 and IRAS F13469-H5833, 
for which the Ha data have low S/N, Mrk23 1 which is severely 
contaminated by an AGN, and IRAS F13342-H3932, which has 
another AGN and it is observed near the edge of the FoV and 
therefore a reliable size determination cannot be made. In gen- 
eral, the INTEGRAL observations were carried out with bun- 
dle SB2 (i.e. optical fibers of 0.9 arcsec in diameter and FoV 
of 12.3 X 16 arcsec^), except for a few cases (see Table 2). The 
spectral resolution is about 6 A. Details about the observations, 
reductions, and calibrations of these data can be found in Garcia- 
Maiin et al (2009a, hereafter GM09). 

The VIMOS sample is presented in Rodriguez-Zauiin et 
al (2011) and is drawn from the IRAS Revised Bright Galaxy 
Sample (Sanders et al. 2003), the IRAS 1 Jy sample of ULIRGs 
(Kim et al. 1998), and the HST/WFPC2 snapshot sample of 
bright ULIRGs (ID 6346 PI: K.Bome). These are mainly south- 
ern objects. For the present analysis, we exclude F08424-3130 
and F12596-1529 because they are not well-covered by the in- 
strument FoV. We used the high resolution mode with the HR- 
orange grating, which provides a spectral resolution of about 2 
A. After combining four dither pointings, the total FoV is about 
30 X 30 arcsec^, with a spaxel scale of 0.67 arcsec (square). 
Details about the observations, data reduction, and calibration 
can be found in Monreal-Ibero et al. (2010) and Rodriguez- 
Zaurin et al. (201 1, hereafter RZl 1). 

In summary, the whole sample consists of 54 systems (66 
individual galaxiesfl 32 LIRGs, and 22 ULIRGs. Their mean 
(median) distance is 226 (137) Mpc, covering a range from 40.4 
Mpc ( z= 0.0093) to 898 Mpc (z= 0.185). Their infrared lumi- 
nosity spreads over the range IO'^'^Lq < L,v < IO'^'^Lq. It also 
includes all types of nuclear excitations and interaction phases, 
and therefore it should be representative of the general properties 
of local U/LIRGs. 

3. Hot emitting region extension in local U/LIRGs 

Leaving aside AGN and shock effects, the Ha emission traces 
the (moderately obscured) gaseous regions ionized by young 
massive stars. Obtaining Ha emission maps is generally costly 
as it requires the use of some kind of 3D observational technique. 
For the present study, we obtain the Ha maps from the IFS data 
cubes (see Sec. 2) after fitting the Ha- [Nil] complex in the indi- 
vidual spectra associated with each spaxel (GM09, RZl 1). From 
these maps, we derive half-light radii (i.e. ri/2). 

A method commonly used to infer half-light radii is based on 
fitting the observed flux distribution to a galaxy model assuming 
some standard surface brightness profiles (e.g. GALFIT, Peng et 
al. 2010). This method is accurate as long as the model is a good 
representation of the actual galaxy flux distribution. This is usu- 
ally the case when studying the starlight of the stellar population 
traced by the rest-frame optical and near-infrared continuum in 
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normal galaxies. However, the star forming regions traced by 
their Ha emission have usually irregular, clumpy, ring-shaped, 
or other peculiar morphologies that cannot be adequately mod- 
eled with standard profiles. This is especially true for interacting 
and merger systems such as U/LIRGs observed at relatively high 
spatial resolutions (i.e. sub-kpc scales). In these cases, the half- 
light radii can be obtained from the curve-of-growth (CoG) of 
the flux in increasingly larger apertures. We follow this approach 
for the local sample (Table 2), as most previously studied mas- 
sive high-z galaxies have determinations based on this method. 
In addition, this method depends less on angular resolution ef- 
fects than others (Forster Schreiber et al, 2009). However, it is 
important to realize that the ri/2 (CoG) of the Ha emitting re- 
gion is not necessarily a measure of the absolute size of the star- 
burst (if such a concept exists for a clumpy flux-uneven struc- 
ture). These measurements should only be considered within the 
constraints imposed by the tracer, angular resolution, and radius 
definition (e.g. CoG) used. These Ha determinations for local 
samples are of relevance because similar measurements are be- 
ing obtained for high-z samples, hence a direct comparison is 
possible. 

In Table 2, we also include size estimates obtained using an 
alternative method (which we refer to as the A/2 method) that 
computes ri/2= VA/tt, where A is the angular extent of the min- 
imum number of pixels (or spaxels) encompassing half of the 
flux. The results of this method do not depend on the location of 
individual emitting regions across the Fo V, but on their actual ex- 
tent. The method does not require knowledge of the galaxy cen- 
ter, which in some cases is off'set between the different galaxy 
components (i.e. Ha and the continuum emission; e.g. Garcia- 
Marin et al. 2009b). This method has been previously used in 
high-z samples (e.g. Erb et al. 2004), though the CoG is more 
commonly applied. The CoG and A/2 methods provide similar 
results for a smooth radially decaying flux emission, but differ 
for highly structured objects. The detailed comparison among 
different tracers and methods for estimating sizes in SFGs will 
be made somewhere else. Here we derive Ha sizes based on the 
CoG method in order to homogenize the comparison with high-z 
studies. 

The sizes obtained directly from the Ha images were trans- 
formed into intrinsic (i.e. deconvolved) sizes by subtracting in 
quadrature the PSF (Table 2), as done for high-z samples. In 
Figure 1 , we present these results as a function of the infrared 
luminosity. We distinguish between interacting and pre- coa- 
lescence merger systems (triangles), post-coalescence mergers 
(squares), and isolated disk galaxies (circles). Five objects have 
half-light radii equal to or smaller than the PSF, and are rep- 
resented as upper limits. Measured sizes smaller than the PSF 
can be obtained as a consequence of seeing fluctuations. For a 
relatively large fraction of LIRGs only a lower limit to the size 
could be obtained. These are objects that, after visual inspec- 
tion of their Ha maps, have clear evidence that a fraction of the 
emission is lying outside the IFU frame. 

The sample shows a wide range of sizes, from ri/2 ~ 0.2 up 
to 7 kpc (see Tables 1 and 2). To analyze the size distribution, we 
restrict the sample to objects with log (L,>/Lo) > 11 .4 because the 
problems due to the limited FoV are significantly smaller This 
group contains a large fraction of objects with small sizes: 39 
% (16/41) have r 1/2 < 1 kpc, and 68 % (28/41) ri/2 < 2 kpc. 
However, a sizable fraction (13/41) of objects have large Ha 
emitting regions (i.e. ri/2 > 2 kpc). Most of these large objects 
are ULIRGs in a pre-coalescense merger phase, with nuclear 
separations ranging from 1.5 kpc to 14 kpc (GM09a). The only 
exception is the ULIRG IRAS 1 1087+5351, which has a nuclear 



separation of 1.5 kpc, at the border considered by GM09a to 
distinguish pre- and post-coalescence systems. The characteri- 
zation of the merger phase in these ULIRGs is secure in most of 
the cases thanks to HST imaging, which also reveals large (30-50 
kpc) envelopes associated with the old stellar component and/or 
with prominent tidal structures (GM09a; see also Sec. 4. 1). 

For many objects with log (Lu/L©) < 1 1 .4 we could only in- 
fer lower limits to the sizes of the Ha emitting region because 
of the limited FoV of the IFS instruments. This luminosity range 
includes objects with a wide range of morphologies (e.g. Arribas 
et al. 2004), and extended Ha emission is not necessarily associ- 
ated with mergers, but with isolated disks as well. We note that, 
on the basis of NICMOS imaging of a distance limited sample 
of 30 LIRGs, Alonso-Hen-ero et al. (2006) found that about half 
of the sample have compact (1-2 kpc) Paa emission with a high 
surface brightness, while for the remaining half the emission ex- 
tends over scales of 3-7 kpc and larger 
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Fig. 1. Intrinsic Ha half-light radii (CoG, method) of local 
U/LIRGs as a function of the infrared luminosity. The symbols 
indicate different morphological classes: triangles are interact- 
ing and pre-coalescence systems, squares are post- coalescense 
mergers, and circles are isolated disk galaxies. Green symbols 
are objects included in the subsample for comparison with high- 
z. The individual galaxies of the system IRAS 06259-4708 (C 
and S), which also belong to this subsample, are not included in 
the plot since no values for their L,> exist. However, consider- 
ing the total luminosity of the system (log(L,>/LQ)= 11.91) and 
their Ha luminosities, they are likely LIRGs. Open symbols are 
used for LIRGs, and solid ones for ULIRGs. Values encircled 
in blue indicate targets with evidence of hosting a (weak) AGN 
according to their optical spectra (see Table 2). For unresolved 
systems, upper limits are represented. For some low luminos- 
ity objects, the values represented are lower limits owing to the 
limited FoV. This is illustrated by the arrows, which indicate a 
factor of 1.5 increase in size. The blue cross in the lower-right 
corner indicates the typical errors of an object with ri/2 ~ 2 kpc 
(for detailed individual errors see Table 2). See text. 
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Table 1. Ha half-light radius and luminosity and SFR surface density for different local U/LIRGs and high-z samples. 



Sample 


n 




ri/2(Ha) 
(kpc) 


(er 


■g s^'kpc"-) 


(Mo 
Median 


^SFR 

yr^'kpc"^) 
Range 


Comments 


Local U/LlRGs 


Median 


Range 


Median 


Range 


i 


log(L,>/Lo)> 11.4 (all) 
log(L,>/Lo)>11.4(noAGN)) 
log(U/Lo)< 11.4 (all) 


41 
33 
20 


1.1 
1.0 

>1.7 


0.2-7.1 
0.2-7.1 
0.2 - >3.7 


41.6 
41.4 
<40.6 


40.5-43.1 
40.5-43.1 
<39.6-42.4 


3.1 

2.2 
<0.3 


0.2-91. 
0.2-91. 
<.03-18. 


ii 
iii 

iv 


Local U/LlRGs subsample 
for comparison with high-z 
















V 


All 

No AGNs 

No pre-coalescence 

Only pre-coalescence 


26 
20 
10 
16 


1.9 
1.6 
1.0 
2.9 


0.3-7.1 
<0.5-7.1 
0.3- 1.9 
0.4-7.1 


41.6 
41.7 

42.2 
41.2 


40.5-43.1 
40.5-43.1 
41.1-43.1 
40.5-42.6 


3.1 
3.9 
14. 
1.1 


0.2-91. 
0.2-91. 
1.1-91. 

0.2-34. 


iii 
vi 
vii 


High-z SGFs 
















VUl 


All 

IR selected 


81 

34 


2.8 
3.4 


0.6-7.5 
1.3-7.5 


41.5 
41.5 


40.2-42.8 
40.2-42.6 


2.7 
2.7 


0.1-50. 
0.1-34. 


ix 



Notes. Main columns are: (1) Sample: For the characteristics of the different samples considered see comments in column (6) and main text; (2) 
n: number of objects in each sample; (3) ri/2(HQ-): intrinsic (deconvolved) Ha half-light radius; (4) logarithm of the dereddened Ha luminosity 
surface density within the half-light radius, ri/2(HQ'). Note that Z/y„ = L(HQ')/(2x;rr^,2), where LCHa) is the total Ha luminosity and the factor 2 
takes into account that within the half-light radius only half of the flux is included; (5) Star formation rate surface density within ri/2(Ha), obtained 
from S„„ following Kennicutt (1998); (6) Comments where: (i) Sample observed via IFS with INTEGRAL and VlMOS-lFUs. Mrk 273 is not 
included in the luminosity calculations, as no IFS calibrated data exists (GM09b). For Arp 220, the Ha luminosity was obtained from Colina et 
al. (2004), while for the rest of the sample they come from GM09b and RZl 1 (see text). This sample is divided at log(L„-/Lo)=l 1.4 because ri/2 
values for the subsample above this limit are virtually unaffected by FoV, while for a large fraction of objects below this luminosity only lower 
limits to their size could be obtained (see text, and Fig. 1). (ii) Galaxies C and S of IRAS 06259-4708 are not included as it is uncertain whether 
they belong to this group as no reliable individual values for their infrared luminosity exits. In any case, the reported statistical values do not 
change by their inclusion, (iii) Objects with evidence of hosting an AGN have been removed from this sample, (iv) For a large fraction of objects 
in this luminosity range, the ri/2 determinations are significantly affected by the small FoV of the IFS instruments used and, therefore, these values 
should be only considered as lower limits, (v) Local U/LIRGs with (dereddened) L(Ha) > 10"*^ ergs"' and observed on similar linear scales as 
the selected high-z samples, (vi) This set excludes objects classified as interacting and pre-coalescence systems. IRAS Fl 1087+5355, which is at 
the border of the definition of a pre-coalescence system (see GM09a) is also excluded, (vii) Only interacting and pre-coalescence systems (IRAS 
Fl 1087-1-5355 is also included), (viii) High-z SFGs observed in Ha with IFS systems with high angular resolutions. They include the works of: 
Forster Schreiber et al. (2009), Law et al. (2009), Wright et al. (2009), and Wisnioski et al. (201 1). (ix) This sample includes only the IR-selected 
sources of the SINS sample of Forster Schreiber et al. (2009). 



3.1. Reddening effects 

A radial variation in the extinction modifies the flux distri- 
bution and therefore affects the radius determinations. The 
U/LIRGs have large amounts of dust in their nuclear regions 
(e.g. GM09b), and therefore extinction-corrected half-light radii 
should be smaller than uncorrected ones. However, detailed red- 
dening corrections of the radius determinations are beyond our 
current possibilities. Aside from the Balmer decrement provid- 
ing only a partial estimate of the actual extinction (e.g. Alonso- 
Herrero et al. 2006), the 2D reddening structure is mostly un- 
known. In general, H/3, and therefore the Balmer decrement, can 
only be obtained for a few objects observed with INTEGRAL 
(GM09a), and is restricted to the innermost regions where H/3 is 
detected with sufficient S/N. Hence, this index cannot be used to 
map the extinction over a sufficiently large area to correct indi- 
vidually the radii. However, these IFS data provide reliable in- 
dividual nuclear extinction values, as well as some general con- 
straints on the extranuclear region. In particular, GM09b found 
that extranuclear extinction for ULIRGs is patchy, with a large 
scatter among objects and, in general, undergoes a radial decay. 
To obtain a coarse estimate of the importance of the extinc- 
tion effects on our determinations, we correct our images with a 
simple model of extinction consisting of a linear decay from the 
nuclear value of Av (inferred from the Balmer decrement) up 
to Av=0 in the outermost regions of the extended Ha emission. 
This implies gradients in the range of 0.15 - 1 visual magnitudes 
of extinction per kpc. This model has two main limitations. First, 



it does not capture the patchy nature of the extinction. Second, 
it only considers the extinction inferred from the Balmer decre- 
ment, which as mentioned above gives only a lower limit to the 
true obscuration. In principle, the limitations of using the Balmer 
decrement should be greater for objects with high extinction in 
their nuclei. However, this model can provide us with a refer- 
ence for evaluating the importance of extinction effects. When 
applying this correction to the sample 0, we found that ri/2(Av- 
corrected) / ri/2 = 0.72 + 0. 15 (also see Fig. 2 for the distribution 
of values), and therefore a mean reduction in size of about 25- 
30 percent. We attempted to analyze a possible correlation be- 
tween these correction factors and L,>, but the low statistics and 
the large scatter among the objects prevent us from reaching any 
firm conclusion. 

In this context, it is also relevant to mention that Alonso- 
Herrero et al. (2009) found that ground-based IFS Ha and 
NICMOS Paa emissions of LIRGs have similar morphologies, 
suggesting that the extinction effects on Ha are not severe, ex- 
cept in the very nuclear regions. 

As for the continuum emission, Veilleux et al. (2002) ob- 
tained ground-based optical R and near-infrared K imaging for 
most of the 118 ULIRGs from the IRAS 1 Jy sample, and in- 
ferred mean half-light radii of 4.80 + 1.37 and 3.48 ± 1.39 kpc, 
respectively. This difference (38%) cannot be attributed totally to 
the difference in extinction between the two bands, since looking 
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0.5 1 

ri/2(Av-corr)/ri/2 

Fig. 2. Distribution of the ratio between the half-Hght radii ob- 
tained after correcting from extinction with a simple model (see 
text), and the ones derived directly without correction. 



at the original images (Kim et al. 2002) it is clear that the K-band 
images are considerably shallower, which make their size deter- 
minations relatively smaller 

Summarizing, though a detailed size correction for redden- 
ing is impossible, our simple model as well as other empirical 
results (Alonso-Herrero et al. 2006 and Veilleux et al. 2002) do 
not suggest that the (observed) reddening has a severe impact on 
the radius determinations. For the remainder part of the paper, 
we use uncorrected size measurements, and therefore they refer 
to the extension of the unobscured Ha emitting region. 

3.2. Ha versus MIR emission extension 

A detailed comparison of our Ha sizes with previous MIR mea- 
surements is difficult because of the different methodology, an- 
gular resolution, and size definition used by the different works, 
and is beyond the goals of this paper However, the general con- 
clusions show relatively good agreement. On the one hand, Diaz- 
Santos et al. (2010), who used long-slit IRS spectra to estimate 
the fraction of extended emission, showed that the MIR contin- 
uum (i.e. 13.2 yum) in LIRGs originates on scales of up to lOkpc, 
with a mean size of the cores of 2.6 kpc. These figures are con- 
sistent with our findings, especially taking into account the dif- 
ferent size definition and angular resolution. They also find that 
for ULIRGs the MIR emission is more compact. Our data do not 
allow us to study in detail the dependence of size on luminosity, 
because for many of the low luminosity objects we could only 
infer lower limits to their size. However, the comparison of the 
median values of our high and low luminosity bins (see Table 1 ), 
suggests a similar behavior to that of Diaz-Santos et al. (2010). 
Finally, they also found that the compactness strongly increases 
in objects classified as mergers in their final stage of interaction. 
This is also clearly observed in our sample when comparing pre- 
and post- coalescence systems (Fig. 1). Diaz-Santos et al. (201 1) 
found that the [Nell] 12.8 /im emission is as compact as the con- 
tinuum dust emission, though the polycyclic aromatic hydrocar- 
bon (PAH) emission is more extended. 



In addition, there is evidence that the Ha, Paa, and MIR 
emissions in LIRGs are strongly correlated with each other On 
the one hand, Alonso-Herrero et al. (2006b) and Diaz-Santos et 
al. (2008) showed that the overall morphologies of the MIR 8/im 
emission (produced by thermal continuum from hot dust and by 
a PAH feature) and the Paa emission line of LIRGs are similar. 
Moreover, Alonso-Herrero et al. ( 2009) found similar IFS Ha 
and NICMOS Paa morphologies, both tracing the nuclear emis- 
sion as well as the emission from bright high surface-brightness 
HII regions. In addition, Alonso-Herrero et al. (2006b) showed 
that the fraction of total emission contained in the relatively 
small NICMOS FoV (19"xl9") is similar for Ha and 24yL(m for 
three LIRGs (see their table 7), suggesting further that a large 
difference among the global extension derived from Ha and the 
MIR emission is not expected. 

3.3. AGN effects 

The presence of a bright AGN could affect the size determina- 
tion, as the extra flux in the nuclear regions associated to the 
AGN would reduce the derived half-light radius. From the origi- 
nal INTEGRAL and VIMOS samples, several objects exhibiting 
a strong influence by the AGN on their optical spectra were ex- 
cluded (e.g. Mrk231, IRAS Fl 3342+3932). However, we retain 
other objects with hints of a weak or modest AGN contamination 
(see Table 2), for which the radius measurements should not in 
principle be severely affected. We note that removing all objects 
with evidence of an AGN may also introduce some systematic 
biases in U/LIRG samples because the presence of an AGN cor- 
relates with other object properties (e.g. Veilleux et al. 1999). 
However, we mark these objects in the figures and discuss the 
possible influence of the AGN on specific results. This approach 
is similar to that of other studies of high-z samples (e.g. FS09), 
which are known to be contaminated by AGNs (see also Shapiro 
et al, 2009). 

Fig. 1 shows that those galaxies with evidence of hosting an 
AGN do not have a significantly distinct behavior from the rest 
of the sample. As can also be seen in Table 1, this is particu- 
larly true for the subsample of galaxies with log(L,>/L0) > 1 1 .4, 
and for that used for high-z comparison (see Sec. 4. 1). However, 
for the low luminosity bin (i.e. L/j./Lq) < 11.4), objects with an 
AGN seem to be on average smaller than those without traces of 
activity. This result should be interprted with caution as it may 
be a consequence of the small number statistics. 

In this context, we consider the results of Alonso-Herrero 
et al (2012) who, based on Spitzer data, found that only 8% 
of local LIRGs have a significant AGN bolometric contribution 
Lboi[AGN]fLir > 0.25. For ULIRGs, the effects are expected to 
be somewhat larger as a consequence of the well-known trend 
of increasing AGN significance with bolometric luminosity (e.g. 
Veilleux et al. 1999, Nardini et al. 2010). 

4. Direct comparison with high-z samples 

4.1. IFS based Ha emission size of local and high-z SFGs 
samples 

One should ideally use not only the same tracer, but also data 
of similar linear resolutions when comparing the sizes of SF 
regions of different galaxy samples. This is in general difficult 
when comparing local and high-z samples owing to the differ- 
ent observed-frame wavelengths of the tracer, and the generally 
different linear resolutions on target. 
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Major efforts have been made to characterize high-z SFG 
galaxies by different groups, and a number of samples of up 
to z ~ 2.5 have been observed in Ha (e.g. Erb et al. 2006a,b, 
Genzel et al. 2008, Wright et al. 2009, Law et al. 2009, Forster- 
Schreiber et al. 2009, 2011, Epinat et al. 2009, Jones et al. 2010, 
Wisnioski et al. 201 1, Nelson et al. 2012, Sobral et al. 2012 and 
references therein). Some of these studies are based on near-IR 
IPS with high angular resolution (i.e. either AO-assisted or under 
good seeing conditions), providing sub- to few kpc linear resolu- 
tions. Our seeing-limited optical IPS observations probe the Ha 
emission of local U/LIRGs at linear resolutions ranging approx- 
imately from 0.2 kpc to 4 kpc. Therefore, there is an overlapping 
fraction of local and distant galaxies for which Ha sizes can be 
compared in a rather homogeneous way, with of similar linear 
resolutions. 

To select from our sample objects well-suited to a proper 
comparison, we imposed the following two conditions: that i) 
their linear resolution, and ii) their (dereddened) Ha luminos- 
ity should be within the range of the high-z samples. Regarding 
the spatial resolution, we restricted the local sample to galax- 
ies observed on spatial scales larger than 400 pc / pixel. This 
physical scale is equivalent to that obtained with a 50 mas 
/pixel at z ~ 2 (i.e. the finest plate scale used so far). The sec- 
ond criterion restricts the local comparison sample to galaxies 
with L(Hq') > 10"*^ ergs"'. After imposing these two conditions, 
we ended up with a comparison sample of 26 local systems. 
Eighteen objects are ULIRGs, and six high-luminosity LIRGs 
(i.e. log(L,>/L0)>11.5). [j Therefore, the comparison sample is 
formed predominantly by our most luminous systems, with a 
majority of ULIRGs. Six objects have evidence of hosting a 
weak AGN based on of their optical spectra (see Table 1)0 

As for the high-z samples, we selected only sources observed 
with IPS in Ha under good angular resolution conditions (i.e. 
with an average PWHM (PSP) of < 0.6 arcsec - 5 kpc at z ~ 2). 
These include the SINS sample (Porster Schreiber et al. 2009, 
- hereafter PS09-, and references therein), and those observed 
by Law et al. (2009), Wright et al. (2009), and Wisnioski et al. 
(201 1). The SINS sample consists of 62 objects at z ~ 2 detected 
in Ha with SINPONI (Eisenhauer et al. 2003) at the VLT. Most 
of the sample consists of infrared-selected objects covering the 
mass range 2x10'' Mq < M* < 3.2x10" Mq, and it is consid- 
ered to be a good representation of massive (M* > 10'° Mq) 
actively star-forming galaxies at that redshift (PS09). Several 
sources were observed with the AO system. Wright et al. (2009) 
observed 6 SPGs at 1.5 < z < 1.7 at the Keck II telescope us- 
ing the near-infrared integral field spectrograph OSIRIS (Larkin 
et al. 2006) on the LGS-AO 0.1" lenslet scale. The objects were 
selected from the rest-frame UV color-selected catalog of Steidel 
et al. (2004), and cover the stellar mass range 2x10^ - 1.6x10'" 
Mq. Using the same instrument Law et al. (2009) detected Ha 
emission in 1 1 SPGs at 2.00 < z < 2.42 from a larger observing 
sample. With stellar masses in the range 10^ - 8x10'" Mq, these 
galaxies are considered to be representative of the mean stellar 



^ The remaining two are the galaxies C and S in the interacting sys- 
tem IRAS F06259-4708, for which no accurate individual L,> values ex- 
ist. However, taking into account the infrared luminosity of the whole 
system log(L„VL0 = 11.91), and their Ha luminosities, they are likely 
LIRGs. 

'' IRAS F07027-6011N, which meets the luminosity and resolution 
criteria, was not included in this sample as it has an Ho- luminosity that 
is anomalously high (2.92 xlO'*~ ergs"') for its relatively low infrared 
luminosity (log (L„/Lo)= 1 1.02). This, together with a broad Ho- emis- 
sion line and a high [NII]/Ha ratio, is a clear indication of significant 
AGN contamination. 



mass of star-forming galaxies at similar redshifts. The selected 
spaxel scale was 50 mas but to increase the S/N, the data-cubes 
were smoothed with a Gaussian of a typical PWHM of 150 mas. 
Using the same instrument, Wisnioski et al. 2011 obtained AO- 
assisted IPS data for a sample of 13 SPGs from the UV-selected 
WiggleZ Dark Energy Survey, with strong [Oil] emission lines 
and 1.28 <z < 1.46. Their stellar masses are in the range 6.3x10^ 
- 5x10"Mo, so are at the high end of the stellar mass distribu- 
tion probed by IPS samples. As a whole, these samples should 
therefore be representative of the intermediate - to-massive SPGs 
at z ~ 2. 

In Figure 3, we compare the intrinsic (deconvolved) Ha sizes 
obtained for the local U/LIRGs with those for the high-z samples 
from the Uterature. We note that this is a rather homogeneous and 
direct comparison as it is done using the same tracer (i.e. Ha), 
the same technique (i.e. IPS), and similar linear resolutions (i.e. 
~ one to few kpc scale). In addition, we used the same method 
to infer half-light radii as that used for most of high-z determina- 
tions (i.e. Curve-of-Growth). In the bottom panel, the extinction- 
corrected L(Ha)|3is shown as a function of ri/2 (Ha). To discuss 
potential biases caused by the limited angular resolution, we plot 
in the top panel of the figure the radii as a function of the spatial 
scale of the resolution element used in the different observations. 
This corresponds to the linear coverage on target of one spaxel or 
half the PWHM of the PSP, whichever dominates the actual res- 
olution. 13 We note that this is merely an indication of the actual 
spatial resolution, which also depends on the changing seeing 
conditions. In the top panel, we also plot the region out of reach 
owing to the limited angular resolution (i.e. left of the dashed 
line).a 



' Despite the interest of SMGs galaxies to the study U/LIRGs (e.g. 
Tacconi et al. 2006, 2008; Swinbank et al. 2004), there are very few IPS 
Ha data sets for these objects (though see Tecza et al. 2004; Nesvadba 
et al. 2007; Swinbank et al. 2006). Unfortunately, the AGN contamina- 
tion and the lack of Ha half-light radius determinations in most cases, 
prevent us from making a direct and homogeneous size comparison for 
this class. 

* Por the high-z samples, Av values are obtained from the literature 
SPD fittings to broad-band global magnitudes (typically, UV, optical 
and near-IR, and, in some cases, IRAC and MIPS data. Por further 
details, see references below). These Av(SPD) were transformed to 
Av(nebular) following Calzetti et al. (2000), who measured P„j,;;n,-(B- 
V) = (0.44 + 0,03)P„rf„to,(B-V) (i.e. Av(nebular)=Av(SPD)/0.44). 
Therefore, for those works that used Av(SPD) to correct the observed 
luminosities (Wright et al. 2009, Law et al. 2009, Wisnioski et al 201 1), 
this correction was recalculated by considering Av(nebular). Por the lo- 
cal sample, Av values are mainly derived from the Balmer decrement 
within apertures of about 1-3 kpc (see GM09, RZll). Por most of the 
VIMOS sample, the Ha luminosities were obtained directly from the 
fluxes given by RZll. Por the objects that were not corrected for red- 
dening in RZll owing to the lack of measurements of their Balmer 
decrement, we took P(B-V)=0.9. This is the mean value obtained for 
the rest of the sample, and agrees well with the one obtained by Veilleux 
et al. (1999) for their sample of ULIRGs (using an extraction aperture 
of 2 kpc). Por the INTEGRAL sample, we took the extinction values 
from GM09, which were calculated within typical apertures of 2-3kpc. 
That both, Av(nebular) and Av(SED) values are based on flux-weighted 
measurements reduces the efl^ects of the different physical scales used. 

' Por our local IPS data we took the spaxel size, for Wright et al 2009 
data we consider 0.1 arcsec, for Law et al. 2009 data we consider half 
the PWHM of the PSP after smoothing (see their Table 1), for Wisnioski 
et al. 201 1 data 0.05 arcsec, and for SINS data half of the PSP PWHM, 
and for objects without seeing measurements a mean value of 0.275 
arcsec. 

** Due to its approximative nature some values can appear in the for- 
bidden region. 
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Fig. 3. Deconvolved half-light radius derived using the Curve of Growth method as a function of the reddening corrected Ha 
luminosity (lower panel) and the spatial sampling (upper panel). For the local U/LIRG sample (green symbols), we used the same 
code as in Fig. 1 . Only the local U/LIRGs with (reddening corrected) Ha luminosities above 10^^ erg s ' and observed with a spatial 
resolution similar to that achieved for distant galaxies (i.e. sampling > 400 pc spx"') are considered. The high-z SFG samples are 
represented with the following symbol code: Red circles, SINS sample (Forster-Schreiber et al. 2009); magenta crosses. Law et al 
(2009); red triangles, Wright et al. (2009); orange circles, Wisnioski et al. (2011). Typical errors in ri/2 range from 30 to 50 percent 
(for detailed individual errors see Table 2). The dashed line in the upper panel approximately defines the region out of reach owing 
to the limited resolution (i.e. left of the line). The cases for which the derived radii were equal to or smaller than those of a point 
source are shown as upper limits. The histograms for the local (green) and high-z (red) samples are normalized to the total number 
of objects in each sample. 



From the comparison of distant massive SFGs and local 
U/LIRG in Figure 3, we can draw several results (also see Table 

1). 

First, local U/LIRGs have a similar range of sizes as high-z 
galaxies. Some local U/LIRGs behave in this plot as the small 
compact objects observed by Law et al (2009), while the rest 
covers the area defined by the other high-z samples. 

Second, according to the size histograms (bottom panel), 
small objects (i.e. ri/2 < 1 kpc) exist with a higher frequency 



in the local sample of U/LIRGsQ In principle, one could think 
that this is a selection effect since such compact objects cannot 
be resolved with SINFONI under seeing-limited conditions, and 
this is the observing set-up for the majority of the high-z galaxies 
in the plot. As indicated by the dashed line in the top panel of the 
figure, small distant galaxies cannot be clearly probed in typical 



' A Kolmogorov-SmimofF test indicates that the two size distribu- 
tions do not match with a probability of 99 % or more. 



Please give a shorter version witli: \authorrunning and/or \titilerunning prior to \maketitle 



seeing resolutions. Interestingly, while we find four unresolved 
objects in the local sample, all the high-z galaxies can be re- 
solved. [3 All this supports the result that galaxies with compact 
Ha emitting regions, such as those found locally, are less fre- 
quent at high-z. The distributions indeed show that while at high- 
z fewer than 3 percent (2 out of 81) have ri/2< 1 kpc, and fewer 
than 25 percent (20 / 81) have ri/2< 2 kpc, much higher percent- 
ages of 23 % (6/26) and 58 % (15/26), respectively, are found for 
the local sample. We note that the difference between the local 
and distant galaxy size distributions would be even clearer if the 
results of Law et al. 2009 were biased towards small radii as a 
consequence of the low sensitivity of the instrumental configura- 
tion (see Table 1). This possibility has been suggested (Law et al. 
2009, FS09), but it is still unclear whether the sizes found could 
be the result of the intrinsic properties of these galaxies, which, 
on average, are somewhat less massive than the other samples. 

We now check for possible systematic effects introduced by 
AGNs (see Sec. 3.3). In Fig. 2, we observe that the systems 
with evidence of hosting an (optically weak) AGNQj do not af- 
fect the general size distribution (also see Table 1). High-z sam- 
ples also contain known AGNs (e.g. SINS, FS09; Shapiro et al. 
2009). However, the five AGNs in the SINS sample (Q1623- 
BX663, K20-ID5, D3a-7144, D3a-15504, ZC-1 101592) have 
normal sizes with ri/2(HcK) ranging from 3.4 kpc to 5.0 kpc (see 
their Table 6), hence the effects of the AGNs in this comparison 
are estimated to be small. 

Another interesting result drawn from the figure is that for 
a sizable fraction of local U/LlRGs, Ha emitting regions with 
large sizes, similar to those of many high-z objects, do exist. 
In fact, 30 % (8 out of 26) of the objects in the present lo- 
cal sample have half-light radii larger than 3 kpc. All of these 
systems are ULIRGs in a pre-coalescence merger phase, as in- 
dicated in Section 3.1. None is a single nucleus object, either 
a disk or a post-merger However, we note that these objects 
were classified as pre-coalescence mergers thanks to relatively 
deep broad-band continuum imaging obtained from the ground 
and, in most of the cases, with HST as well. In general, there 
is no imaging of similar quality for the high-z samples, hence 
their structure and dynamical phase are not known with a sim- 
ilar level of detail. However, Forster-Schreiber et al. (201 1) ob- 
tained deep rest-frame (continuum) optical high-angular resolu- 
tion HST NICMOS imaging for 6 SINS objects, and compared 
it with the Ha emission from their SINFONI IFS data. One of 
these objects (BX528) had been previously classified as a merger 
based on its kinematics (FS09), and has a double-nucleus irreg- 
ular continuum emission (see their Fig. 1), with an associated 
total radius of 4.86 kpc (or 3.18 kpc and 3.57 kpc for the indi- 
vidual components). This object has an ri/2 (Ha) of 4.6 kpc, so it 
has the typical characteristics of one of the pre-coalescence local 
U/LIRGs. The structure of the Ha emission for high-z samples 
(e.g. FS09) often shows a clumpy and irregular structure similar 
to that found in the Ha maps of local pre-coalescence U/LIRGs. 
Furthermore, on the basis of kinematic arguments FS09 found 
that 33 % of the SINS galaxies have evidence of being interact- 
ing and merging systems. Kartaltepe et al. (2011) also found a 
significant fraction of mergers, interacting, and irregular galax- 
ies at z~ 2, especially among ULIRGs. Therefore, it is likely 



that pre-coalescence systems similar to those in the local sam- 
ple also exist in the high-z samples. Interestingly, the median 
properties of the whole high-z sample are similar to those of 
pre-coalescence local U/LIRG systems (Table 1). If the pre- 
coalescence systems are excluded from the local U/LIRGs sam- 
ple, the high-z SF galaxies appear to be larger than the local 
U/LIRGs by a factor of about 3 (see Table 1 for specific values). 
The depth of the observations is another factor that can affect 
the comparison between the local and distant samples. It is diffi- 
cult to estimate the depth of the different observations owing to 
the many factors affecting its calculation, some of which are dif- 
ficult to handle (e.g. seeing fluctuations, wavelength-dependent 
infrared background). However, if the observed integrated Ha 
fluxes of the local U/LIRGs are transformed into surface bright- 
ness and redshifted to z=2.2, there are several cases with val- 
ues below the ~ l.xlO '^ erg s ' found by FS09 for the faintest 
sources of the SINS sample. The objects below this threshold are 
generally large objects of low surface brightness, which would 
have been missed if they were at z ~ 2. Therefore, the relatively 
shallower observations of the high-z samples strengthen the con- 
clusion that local U/LIRGs are on average intrinsically smaller 
than high-z populations. 

4.2. Ha and SF surface densities in star-forming galaxies at 
low and high redshifts 
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Fig. 4. Ha luminosity density as a function of the Ha luminos- 
ity for local U/LIRGs and distant massive SFG with a similar 
luminosity range and observed under similar linear resolution 
conditions. Symbols are the same as in Fig. 3. The histograms 
for the local (green) and high-z (red) samples are normalized to 
the total number of objects in each sample. The blue cross in the 
lower-right comer indicates the typical errors. 



'" Note that SINS objects D3a-7429 and GMASS-1146, for which 
only an upper limit to their sizes were reported by FS09, were observed 
under relatively bad seeing conditions, hence implying a spatial sam- 
pling not probed by the plot. 

" Note that those objects with evidence that a strong AGN affects the 
size determination were removed from the comparison sample. 



Recent works (e.g. lono et al. 2009, Rujopakam et al. 201 1) 
have suggested that the overall size of the star- forming region 
in local U/LIRGs is significantly smaller than in both high-z 
U/LIRGs and SMGs, and local lower-luminosity star-forming 
galaxies. Rujopakarn et al. (201 1) suggested that the star forma- 
tion in local U/LIRGs proceeds with extremely high SF surface 
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densities, in some cases by factors of 1 000 or higher than in other 
SFGs, and therefore that they are driven by a unique process. 
These resuhs are based (1) on the determination of the physi- 
cal sizes using different spectral features for the local and high- 
z star-forming galaxies, and (2) on the assumption that these 
different features give a similar measure of the size of the SF 
region. The available spectral features trace different physical 
mechanisms and components of the ISM. The most commonly 
used tracers are the Ha and Paff lines (warm ionized gas), CO 
lines (cold molecular gas), 24//m continuum emission (hot dust 
emission), and radio continuum (non-thermal synchrotron emis- 
sion). Each of these tracers is affected by its own biases (see 
discussion in Genzel et al. 2010) hence an effort should be made 
to establish the size measurements, and therefore the SF surface 
density, in a homogeneous way when comparing galaxy sam- 
ples. 

To compare our local U/LIRGs with high-z samples, we con- 
verted the extinction-corrected L(Ha) luminosities into Ha lu- 
minosity surface densities within ri/2 (^ho- L(HQ')/27rrj „) us- 
ing the Ha half-light radius derived from the same IFS data 
(Section 3.1). Therefore, both the half-light radius and the lu- 
minosity were inferred from the same set of data, namely the 
IFS-based Ha emission maps. In Figure 4, we make a direct and 
homogeneous comparison between our local U/LIRGs and high- 
z massive SF galaxies covering the same Ha luminosity range, 
and with similar IFS Ha measurements available from the liter- 
ature (see Section 4.1; Table 1). 

Local U/LIRGs and high-z SF galaxies cover a similar range 
in Ha luminosity surface density spanning over at least three or- 
ders of magnitudes from 2//„ of lO'*" to more than 10'*^ erg s"' 
kpc"^. This would correspond to Hsfr between 0.08 and 80 Mq 
yr ' kpc"^ according to Kennicutt (1998). However, while local 
U/LIRGs are more evenly distributed towards high surface den- 
sity values (i.e. "Lna > 10"^^ erg s"' kpc"^), the majority of high-z 
galaxies have moderate surface densities (i.e. £//« between 10'*' 
and 10"^^ erg s"' kpc"^). 

The data in Fig. 4 are also consistent with a trend in the sense 
that the more luminous galaxies also have higher luminosity sur- 
face densities. This linear behavior is more clearly observed in 
the the high-z sample (correlation coefficient = 0.60) than lo- 
cally (0.43). The local galaxies with Ha surface densities that 
are significantly higher than the high-z sample (i.e. > 2cr from 
the mean linear fit) represent about 20% of the sample. If the 
Law et al. (2009) sample (indicated by crosses in Fig. 4) is not 
considered within the high-z pool a stronger correlation is ob- 
tained (0.72), and the percentage of local galaxies that depart 
significantly from the mean high-z behavior increases by up to 
about 40%. 

Only two (optically) identified AGNs are among the local 
U/LIRGs with high Ha surface densities that are therefore sus- 
pected to be contaminated. However, the majority of AGNs are 
in sources with low surface densities, suggesting that they can- 
not be the dominant factor for the local high surface densities. 
Therefore, although the presence of a buried AGN cannot be dis- 
carded (see Section 3.3), the Ha emission is likely due mainly to 
star formation. 

The trend shown in Fig. 4 is equivalent to the one presented 
by Rujopakarn et al. (201 1), based on the total infrared (TIR) lu- 
minosity (i.e. 'Ltir vs. Ltir, see their figure 4). However, while 
in their TIR-based relation most of the the local U/LIRGs largely 
depart from the general trend found in other (local and distant) 
populations, we do not observe this behavior in Ha, except in 
some cases. More specifically, while Rujopakarn et al. (2011) 
found that most of the local U/LIRGs depart by factors of be- 



tween 1 and 3 orders of magnitude in luminosity surface density 
away from other the trend defined by other SFGs samples, our 
Ha based measurements suggests that only a fraction (20-40 %) 
show a clear departure of about one order of magnitude, as can 
be inferred from Fig. 4. 



4.3. Dependence of size on Av and [Nll]/Ha: Comparison of 
locai U/LIRGs witli infrared-seiected SINS galaxies 

As discussed above, the reddening structure may affect the de- 
rived sizes for the Ha emitting region. Therefore, the effects of a 
possible different radial variation in the extinction for local and 
high-z sources should be in principle taken into account when 
comparing both samples. Unfortunately our knowledge of the 
2D dust / reddening structure in ULIRGs is mostly unknown (see 
GM09b, and Sec. 3.1), and we have a complete ignorance of it 
for high-z samples. Therefore a detailed study of this topic is 
well beyond the current possibilities. 

As a first attempt at studying relative effects, we compare 
the Av distributions of the samples under analysis, looking for 
possible correlations with the derived sizes (Fig. 5). No obvious 
correlation size - Av is found for any of the two samples. We note 
that the lack of this correlation does not imply that reddening 
does not affect size, as it is the radial variation of extinction, 
rather than its global value, the main contributing factor. 

Fig. 5 also shows that there is a clear distinction in the be- 
havior of local and distant Av distributions. In particular, while 
the local sample show a large scatter with values in the range 0.5 
- 6 mag, the high-z sample distribution has a lower mean value 
with a peak at about 1.2 mag, hence on average the high-z sam- 
ple is less extincted. However, we see below that a sub-sample of 
high-z targets (i.e. those that are IR-selected) has a distribution 
of Av similar to the local sample. 

For the [NII]/Ha ratio, there is also a clear difference be- 
tween the local and distant samples, as shown in Fig. 6. Ej 
Similarly to Av, there is again no correlation between this ra- 
tio and size. The difference between the two [NII]/Ha distribu- 
tions may be even larger than suggested by the histograms, as 
several values for the high-z sample are only upper limits, as 
indicated by the arrows in the figure. [3 However, it should be 
noted that [NII]/Ha values are available for only a fraction of 



'- For tlie INTEGRAL data, we use the integrated values for the cir- 
cumnuclear region given by Garcia-Marin (2007), while for the VIMOS 
data we consider the integrated spectra following Rosales-Ortega et al. 
(2011). Hence, typical physical scales are a few kpc, which are similar 
to the high-z samples. We note that these are flux-weighted measure- 
ments. 

" If we interpret the [NIIl/Ha ratio as a measure of metallicity via the 
N2 index (e.g. Pettini and Pagel, 2004), we obtain a mean solar value 
of 12 + log (0/H) = 8.70 (AUende-Prieto et al. 2001; Asplund et al. 
2004; Scott et al. 2009) for the present local sample. This agrees within 
the uncertainties with the findings of Rupke et al. (2008), who derive 
metallicities from the nuclear spectra of a local sample of U/LIRGs us- 
ing different calibrators. For the high-z sample as a whole, an upper 
limit of 8.43 is obtained, so at least a factor of two lower metallicity 
with respect to the solar, in agreement with the previous long-slit re- 
sults of Erb et al. (2006b). We should note that, this ratio might also be 
affected by the presence of an AGN and by shock excitation as found in 
local (e.g. Veilleux et al. 1995, 1999; Yuan et al. 2010; Monreal Ibero et 
al. 2006, 2010) and distant (Shapiro et al. 2009) samples. However, ac- 
cording to the diagnostic diagram proposed by Veilleux and Osterbrock 
(1987), the [NIIJ/Hq- values reported in figure 6 suggest that ionization 
is mainly due to stars (i.e. [NIIl/Ha < 0.6), for a wide range of [OIII]/Hq- 
values. 
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Fig. 5. Av (nebular) as a function of Ha radius. Symbol code 
is the same as in previous figures. The histograms for the lo- 
cal and high-z samples are normalized to the total number of 
objects in each sample. For the high-z samples, Av were ob- 
tained from the literature SED fittings to global magnitudes, 
transformed to Av(nebular) following Calzetti et al. (2000) (i.e. 
Av(nebular)=Av(SED)/0.44). For the local sample, Av values 
were derived from the Balmer decrement within apertures of 
about 1-3 kpc (see GM09, RZl 1). 



the SINS sample and these are strongly biased towards the BX 
(i.e. optically-selected ) objects (Forster Schreiber et al. 2006). 

The differences between the mean Av and [NIII/Hq- values 
of the local and distant samples shown above are dominated 
by the contribution of the UV-selected high-z galaxies clustered 
around the low Av values. If we restrict the comparison to the 
infrared-selected galaxies of the SINS sample (see table 1 of 
FS09), the A,, (nebular) distribution is similar to that of the lo- 
cal sample ( <Av,>=3.2 ± 1.5 in both cases), as it can be seen 
in the lower panel of figure 7. For the three infrared-selected 
/ bright SINS objects with [NII]/Ha measurements (D3a-6004, 
Da3- 1 5504, ZC-78294 1 ; FS09, Genzel et al. 2008), a mean value 
of [NII]/Hq' = 0.34 was found, which is close to the value found 
for the local U/LIRGs and significantly higher than those for the 
rest of the high-z galaxies. In Figure 7, we also compare the size 
and luminosity surface density of the local U/LIRGs with those 
of the infrared-selected objects in the SINS sample. Despite the 
similarity in the mean Av and [NII]/Hq' properties, the size dis- 
tributions are different, with a higher frequency of local objects 
of small size (high surface density). Although the number of ob- 
jects is low for a reliable statistics, a Kolmogorov-Smirnoff test 
indicates that the radius distributions do not match with a prob- 
ability of 99% or more. 



Fig. 6. [NIIl/Ho- as a function of Ha radius (see text). Symbol 
code is the same as in previous figures. The histograms for the 
local and high-z samples are normalized to the total number of 
objects in each sample. 
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Fig. 7. Comparison of different properties between the local 
U/LIRGs and the infrared-selected SINS objects from Forster- 
Schreiber et al (2009). The histograms are normalized to the to- 
tal number of values in each sample, which are indicated in the 
upper- right corner of each panel. Note that the metallicity dis- 
tribution for the high-z sample is defined by only three values. 



5. Summary and conclusions 

Using IFS-based Ha emission maps, we have derived the sizes of 
the SF regions for a representative sample of 54 local U/LIRGs 
(66 galaxies). From this initial sample we have selected a sub- 
sample of 26 local U/LIRGs which we have compared with ex- 
isting works of high-z SFGs observed in Ha with near-infrared 



IFS at high angular resolution (i.e. AO-assisted or under good 
seeing conditions) to match the physical size probed locally and 
at high-z. We have also selected these galaxies to have similar 
L(Ha) (i.e. SFRs) as high-z samples. The comparison between 
the basic properties (such as size and luminosity surface den- 
sity) of the galaxies in the local and distant samples has been 
made therefore in a homogeneous way, using the same observ- 
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ing technique (IFS), tracer (Ha), and similar linear resolutions. 
The main conclusions of the present analysis can be summarized 
as follows: 

1- The SF region sizes for the whole local sample of 
U/LIRGs, as derived from Ha, have a wide range of values with 
ri/2(Ha) from 0.2 to 7 kpc. For objects with L,> > lO""* Lg 



■05 



about 2/3 have relatively compact (i.e. ri/2 < 2 kpc) emission. 
For the remaining 1/3 (13/41) with large Ha emitting regions 
(i.e. ri/2 > 2 kpc), 11 are ULIRG in a pre-coalescence merger 
phase with nuclear separations between 1.5 and 14 kpc. For ob- 
jects with L,> < lO'^"* Lo, large Ha emission is not necessarily 
associated to mergers, but to isolated disks as well. 

When we compare the subsample of U/LIRGs with the high- 
z sample (at similar L(Ha) and physical scales) we find the fol- 
lowing additional conclusions: 

2- The Ha size distributions for local U/LIRGs and distant 
massive star forming-galaxies cover a similar range of values. 
However, we have found that there is a higher frequency of com- 
pact objects locally than at high-z. In particular, while objects 
with ri/2 < 2 kpc represent 58 % (15/26) of the low-z subsample, 
they constitute fewer than 25 % (20/81) of the galaxies at high- 
z. The median size of local U/LIRGs is a factor of 1.5 smaller 
than the one for high-z SFGs. This value strongly depends on 
whether only pre-coalescence systems are considered (~ 1) or are 
excluded (~ 3) from the local U/LIRG subsample. These factors 
are on average smaller than the ones reported using a variety of 
other tracers. 

3- Most of the local U/LIRGs and high-z massive SF galax- 
ies cover a similar region in the "Zhu - L( Ha) plane, spanning 
from I.Ha of lO'*" to more than 10"*^ erg s"' kpc"^. About 20-40 
% of the local U/LIRGs studied here show values of I,H,r signif- 
icantly higher (by factors of about 10) than at high-z. These are 
considerably smaller than the much higher factors of ~ 1000 or 
more, recently reported in similar planes (i.e. L(TIR) vs. I.tir)- 

4- When the comparison with the high-z SFGs is restricted 
to the infrared-selected SINS objects, local and distant sam- 
ples have similar properties in terms of global visual extinction 
(and [NII]/Ha). Although with a lower statistical significance, 
the main result of the local vs. high-z comparison remains un- 
changed: the SF region size, on average, is a distinctive factor 
between local U/LIRGs and high-z samples of similar L(Ha) 
(i.e. SFR). 

5- A significant fraction (approximately 1/3) of local 
U/LIRGs, generally pre-coalescence merger systems, are indis- 
tinguishable from typical high-z SFGs galaxies in terms of their 
Ha size and surface brightness. 
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Table 2. Ho- radii for tlie IFS sample of U/LIRGs observed with VIMOS and 
INTEGRAL 



IRAS ID 


Other ID 


Subsample 


ri/z-Ha (A/2) 


Vi/2-na (CoG) 


Comment 


name 






(kpc) 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


VIMOS sample 


F01159-4443N 


ESO-244-G012 


N 


0.35+0.10 


0.36+0.11 


1,2 


F01341-3735N 


ESO-297-G011 


N 


>1.36±0.41 


>1. 89+ 0.57 


2,3 


F01341-3735S 


ESO-297-G012 


N 


0.23+0.11 


0.22+0.11 


1 


F04315-0840 


NGC 1614 


N 


0.51+0.15 


0.51+0.15 




F05 189-2524 




Y 


<0.28 


<0.28 


0,5,6 


F06035-7102 




Y 


3.72+ 1.12 


7.06+2.12 




F06076-2139N 




N 


<0.25 


<0.25 





F06076-2139S 




N 


0.48+ 0.24 


0.45+ 0.22 


1,7 


F06206-6315 




Y 


2.53+ 0.76 


4.43+ 1.33 


5 


F06259-4708N 




N 


0.71+0.21 


0.74+ 0.22 


9 


F06259-4708C 




Y 


1.16+0.35 


1.31+0.39 


8 


F06259-4708S 




Y 


1.45+0.44 


1.93+0.58 


8 


F06295-1735 


ESO-557-G002 


N 


1.91+0.57 


2.46+ 0.74 


2 


F06592-6313 




N 


0.35+0.17 


0.39+ 0.20 


1 


F07027-6011N 


AM 0702-601 


N 


0.47+ 0.24 


0.50+ 0.25 


1,5,6,10 


F07027-6011S 




Y 


0.63+0.19 


0.63+0.19 




F07 160-62 15 


NGC 2369 


N 


>0.54+0.16 


>0.86+0.26 


3 


08355-4944 




N 


0.59+0.18 


0.64+0.19 


2 


F08520-6850 


ESO60-IG016 


Y 


1.02+0.31 


1.04+0.31 


2 


09022-3615 




Y 


1.30+0.39 


1.20+0.36 




F09437+0317N 


IC-563 


N 


>2.20+ 0.66 


>3.22+0.97 


3 


F09437+0317S 


IC-564 


N 


> 1.93+0.58 


>3.08+0.92 


3 


F10015-0614 


NGC-3110 


N 


> 1.95+0.59 


>2.66+0.80 


3 


F10038-3338 


IC2545 


Y 


0.51+0.26 


0.69+ 0.23 


1 


F10257-4339 


NGC 3256 


N 


>1.01+0.30 


>1. 37+ 0.41 


3 


F10409-4556 


ESO-264-G036 


N 


>1.78+0.53 


>3.24+0.97 


3 


F10567-4310 


ESO-264-G057 


N 


>1.61+0.48 


>3.69+l.ll 


2,3 


Fl 1255-4120 


ESO-319-G022 


N 


1.25+0.38 


2.70+0.81 




Fl 1506-3851 


ESO-320-G030 


N 


0.68+ 0.21 


1.05+0.31 




F12043-3140N 


ESO-440-IG 058 


N 


0.24+0.12 


0.18+0.09 


1,7 


F12043-3140S 


ESO-440-IG 058 


N 


>1.38+0.41 


>2.92+0.88 


3 


F121 15-4656 


ESO-267-G030 


N 


1.34+0.40 


1.55+0.46 




12116-5615 




N 


0.28+0.14 


0.29+0.14 


1 


F13001-2339 


ESO-507-G070 


N 


0.88+ 0.26 


0.91+0.27 




F13229-2934 


NGC 5135 


N 


>0.52+0.16 


>0.53+0.16 


3,5 


F14544-4255E 


IC4518 


N 


>0.94+ 0.28 


>1.06+0.32 


3 


F14544-4255W 


IC4518 


N 


0.57+0.17 


0.66+ 0.20 


5 


F17138-1017 




N 


0.56+0.17 


0.63+0.19 




F18093-5744N 




N 


1.10+0.33 


1.56+0.47 




F18093-5744C 




N 


0.35+0.11 


0.34+0.10 


6 


F18093-5744S 




N 


0.84+ 0.25 


1.21+0.36 


7 


F2 1130-4446 




Y 


1.59+0.48 


1.89+0.57 


2 


F21453-3511 


NGC 7130 


N 


> 1.04+ 0.31 


>2.27+0.68 


3,5,6 


F22 132-3705 


IC5179 


N 


> 1.56+ 0.47 


>1.90+0.57 


3 


F22491-1808 




N 


1.70+0.50 


2.30+ 0.70 




F23128-5919 


AM 2312-591 


Y 


2.01+0.60 


3.01+0.90 


5,6 


INTEGRAL sample 


06268+3509 




Y 


3.74+ 1.50 


4.33+ 1.73 




06487+2208 




Y 


0.98+ 0.49 


1.12+0.56 


1,4 


F08572+3515 




Y 


1.85+0.74 


2.70+ 1.08 




Fl 1087+5351 




Y 


2.49+ 1.00 


3.76+ 1.50 


4,6 




Arp299E/IC694'' 


N 


>0.68+ 0.27 


>1.03+0.41 


3 


Fl 1257+5850" 


Arp299W /NGC3690" 


N 


>0.53+0.21 


>0.67+ 0.27 


3,5 


F121 12+0305 




Y 


2.60+ 1.04 


3.48+ 1.39 




F12490-1009 




N 


1.77+0.71 


1.58+0.63 




F13156+0435N 




Y 


1.96+0.78 


1.99+0.80 





Continues on next page 
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IRAS ID 


Other ID 


name 




(1) 


(2) 


F13156+0435S 




F13428+5608 


Mrk 273 


F13536+1836 


Mrk463 


F14060+2919 




F14348-1447 




F15206+3342 




F15250+3609 




F15327+2340 


Arp 220 


F16007+3743 




F17207-00I4 




F18580+6527 





sample 


ri/2-Ha (A/2) 


Ti/2-na (CoG) 


Comment 




(kpc) 


(kpc) 




(3) 


(4) 


(5) 


(6) 


Y 


2.66+ 1.06 


2.68± 1.07 




N 


1.29±0.52 


1.71+0.68 


5 


Y 


<0.44 


<0.44 


0,5,6 


Y 


1.30+0.65 


1.13+0.57 


1 


Y 


3.82+ 1.53 


4.68+ 1.87 




Y 


<1.00 


<1.00 





Y 


<0.48 


<0.48 





N 


0.63+ 0.25 


0.61+0.24 




Y 


3.67+ 1.47 


6.15+2.46 




Y 


0.90+ 0.36 


1.01+0.40 




Y 


2.37+0.95 


2.09+ 0.84 


5 



Notes. Columns: (1) and (2) Identification, (3) Indicates whether or not the object is in the subsample of U/LIRGs with L(H(i') and spatial resolution similar to those of high-z samples 
(see text) (4) Ha radius obtained with the A/2 method, (5) Ho- radius obtained with the Cui've-of-Growth method, (6) comments with the following code: 0:11a half-light radius equal or 
smaller than the PSF (i.e. unresolved), 1: Ho- half-light radius smaller than 1.25 times the radius of the PSF and therefore, with larger uncertainty, 2: Noisy continuum, 3: Limited FoV, 
4: Obsei'ved with INTEGRAL / SBl configuration, 5: Evidence for an AGN according to the emission line ratios (see GM09, and compilation of RZl 1), 6: Evidence for a broad line 
associated to Ha suggesting the presence of an AGN. 7: The L„. for this object is unknown, and therefore it is not included in the analysis. 8: For this object the individual L^,. is not 
known, but taking into account that the whole system has Li, = l 1.91, and (reddening coiTccted) Ha luminosities above 10"^" erg s~' it is likely a LIRG, and it kept for the analysis. 9: For 
this object the individual L„ is not accurate known, but taking into account that the whole system has L„ = 11.91 and according to Spitzer MIPS (24/jm) data it carries a large fraction of 
the IR flux, it is likely a LIRG. 10: This object likely to have a strong AGN contamination so it was not included in the sample for comparison with high-z samples (see text). The values 
reported for the CoG and A/2 methods have been deconvolved in quadrature from the radii associated to stars. We used an averaged value of r|/2 (PSF) = 1.07 arcsec for VIMOS data 
from 7 stars located within the FoV of the science pointings. For the INTEGRAL-SB2 data, that value was scaled according to the spaxel scale, and for SB 1 (which was used on average 
under better seeing conditions) we adopted 0.9 arcsec, as inferred from the Mrk231 obsei'vations. Typical enors are estimated to be 30 percent, except for cases with comments land 2 for 
which are increased to 50 percent. Values limited by FoV (comment 3) should be considered as lower limits. For radii derived from the INTEGRAL data, typical eiTors are estimated to 
be somewhat larger (40 percent), since no direct measurements for the PSF were available. Further details on these galaxies can be found in GM06 (INTEGRAL data) and RZl I (VIMOS 
data). 
" Identification as in Garcia-Man'n et al.(2006). 
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